Interaction ramps in a trapped Bose condensate 



F. E. Zimmer and Masudul Haque 

Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, Germany 

Non-adiabatic interaction ramps are considered for trapped Bose-Einstein condensates. Tlie deviation from 
adiabaticity is characterized through the heating or residual energy produced during the ramp. We find that the 
dependence of the heat on the ramp time is very sensitive to the ramp protocol. We explain features of this 
dependence through a single-parameter effective description based on the dynamics of the condensate size. 
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Adiabaticity is a ubiquitous concept in quantum dynam- 
ics. With the advent of novel non-equilibrium experimental 
possibilities, deviations from adiabaticity in slow parameter 
changes have attracted a lot of attention The ques- 

tion of non-adiabaticity is of fundamental interest, but also has 
practical implications. Many experimental protocols involve 
adiabatically changing a parameter in order to reach a desired 
quantum state. Since non-adiabatic heating can rarely be com- 
pletely avoided, it is essential to understand deviations from 
adiabaticity in slow ramps. In addition, the proposal of adi- 
abatic quantum computation ifioll raises the question of how 
a realistic parameter ramp in a quantum many -particle system 
deviates from adiabaticity. While the effect of quantum criti- 
cal points in the ramp path has been considered in much detail 
IIIII2I], first realizations of such a quantum computer will pre- 
sumably be mesoscopic rather than macroscopic, without true 
quantum critical points. Understanding non-adiabatic ramps 
m finite quantum systems is therefore vital. A few studies of 
ramps in finite and trapped systems have appeared in the very 
recent hterature iSjstl. indicating an emerging recognition of 
the importance of this issue. 

Thus motivated, in this work we consider non-adiabatic ramps 
in the most emblematic physical system in the world of laser- 
cooled atoms, namely, an interacting Bose-Einstein conden- 
sate in a harmonic trap. We consider ramps of the contact 
interaction from an initial value Ui to a final value Uf, occur- 
ring in time scale r. By varying r, we interpolate between 
the limits of the instantaneous quench (r = 0) and the adia- 
batic ramp (t — > 00). For finite r, we study deviations from 
adiabaticity through the heating Q, which we define as the fi- 
nal energy at time t > t minus the ground state enegy of the 
final Hamiltonian. This quantity is also called the residual en- 
ergy or the excess energy , and may be thought of as the 
"friction" due to imperfect adiabaticity jllll . 

We find that the residual energy Q decreases with r as a power 
law Q ^ T^" rather than exponentially, with the exponent ly 
depending on the shape of the ramp. For certain ramp shapes, 
Q{t) has oscillations superposed on top of the power-law de- 
cay. The oscillation frequency is given by the breathing-mode 
frequency set by the harmonic trap. In contrast to most stud- 
ies of the heat function Q{t), we have found a simple phys- 
ical description of the features of Q{t). A natural descrip- 
tion of trapped dynamics is through a variational wave func- 
tion where the extent (radius) of the condensate is treated as 
a time-dependent variational parameter. We find that such a 
"radius dynamics" description reproduces the residual energy 



behavior surprisingly well. 

We use the Gross-Pitaevskii (GP) description lfl3i[l4ll of con- 
densates in an isotropic harmonic trap Vti-{r). The GP energy 
functional is 
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Here r is the radial position variable, and = Jd^r is the 
spatial integral appropriate to the dimensionality D of the sys- 
tem. The time-dependent parameter U is the effective interac- 
tion strength whose relation to the physical interaction is also 
D-dependent (c.f. Ref. II12II for ID). From here on we will use 
trap units, expressing lengths in units of trap oscillator length 
and time in units of inverse trapping frequency. The conden- 
sate dynamics is given by the time-dependent GP equation. 

The GP equation provides an excellent account of many as- 
pects of trapped condensate dynamics. We restrict ourselves 
to GP dynamics, which is already too rich to be studied ex- 
haustively. While physics beyond GP is more important in 
lower dimensions, we will treat ID, 2D and 3D cases on an 
equal footing. The questions we address are conceptually in- 
teresting irrespective of dimensionality. We will also show 
that the dominant effects we encounter and analyze are quite 
insensitive to D. 

Ramp shapes — We analyze ramps of the form 

U{t) = U, + e{t){Uf~Ui)r{tlT). 

The ramp function r{x) starts at r(0) = and ends at 7'(oo) = 
1. We do not require r{x > I) = 1, i.e., the ramps take 
place over time scale r but do not necessarily end at < = r. 
(Contrast, e.g., Ref. fl.) 

Specifically, we consider the following forms for r{x): 



[lA] 
[2A] 
[3A] 
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[IB] 
[2B] 
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Each [A], [B] pair has the same initial behavior, r{x)^x°', but 
the [B] versions have no endpoint kinks. Ramp shapes [lA] 
and [IB] ([2A] and [2B]) are compared in Figure[T^ ([TJ)). 
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FIG. 1: (a&b) Ramp shapes [lA], [IB], & [2A], [2B]. (c-e) Some 
residual energy curves Qij). (c) For ramp shape [lA], Compar- 
ing full Gross-Pitaevskii results (full curves) with single-parameter 
variational results (dotted curves). The three pairs have (Ui,U /) val- 
ues (0,1), (10,100), and (20,0.2), from bottom to top. (d) With 
{Ui,Uf)= (10,100), comparing the discontinuous-derivative ramp 
shapes [lA], [2A], [3A]. (e) With (U^,Uf)= (20,0.2), comparing 
Q(t) for smoothed ramp shapes, [IB], [2B], [3B], from smaller to 
larger slopes. 



Residual energy features — In Figure \T]p-e we present the 
behavior of the heat function Q{t), normaUzed against its 
instantaneouus-quench value Q{t = 0). 
Figures [1]; shows Q(t)/Q(0) for several ramps of type [lA] 
(linear quench). The heating displays oscillations with r, on 
top of a power-law decay. The frequency of these oscillations 
is the same as the frequency of breathing-mode oscillations 
in real time. Also shown, in each case, is Q{t)/Q{0) found 
from a single-parameter variational ansatz where the cloud ra- 
dius is the only variable. The near-perfect agreement indicates 
that the physics of heating in interaction ramps is almost com- 
pletely described by the radius dynamics. In the rest of the 
article, we therefore present results and analysis mostly based 
on the variational description. 

In Figure [TJl, the residual energy curves are compared for the 
ramp shapes r(a;)~2:" with discontinuos derivatives at end- 
points, [1 A], [2A], [3A]. Each curve has an overall power-law 
decay with the same decay exponent, Q{t) ^ r^^. This sug- 
gests that the residual energy for such ramps is primarily set 
by the endpoint kink. Superposed on the power-law decay are 
oscillations, more prominent for the linear quench a = 1 and 
barely visible for a — 3. 



In Figure [T^ we focus on smoothed ramps [IB], [2B], [3B], 
which lead to non-oscillating decay of the residual energy. 
The decay exponent is seen to depend on the power a of 
r{x) ~ x", namely, Q{t) ^ r"^". 

The dimensionality does not affect the decay exponents. The 
Q{t) data shown in Figures[T];-e are for = 1, except for the 
dashed curve in[T^, which is for D = 3 and ramp [3B]. Com- 
parison with the corresponding D ^ 1 solid curve demon- 
strates that the Q(t) behavior is practically identical in differ- 
ent dimensions. 

Single-parameter variational description — We formulate the 
radius description in terms of a Gaussian variational ansatz, 
which for ID is 
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For D > 1 the variational wave function is a product of one 
such Gaussian factor for each dimension. Using this ansatz 
in the GP Lagrangian, we get the evolution equations for the 
variational parameters a{t) and /3{t) lITsll . We could just as 
well use a Thomas-Fermi instead of Gaussian profile; how- 
ever the results are very similar and do not substantially af- 
fect any of the arguments we make in this work. The two 
parameters turn out to be not independent but simply related 
(/3(t) oc dt In cr(t)). There is thus effectively a single dynami- 
cal parameter describing the system, namely the cloud radius 
(T{t). The equation of motion for <t is 
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and the energy is 
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The radius description based on the above two equations 
is suitable for describing breathing-mode oscillations. For 
constant U, trying small-amplitude oscillatory solutions of 
form a = Ra + psm{nt), Eq. ^ yields O ^ V^ZJT2 
for the breathing mode frequency at large U. We also find 
that Rq satisfies the stationary equation for the radius, i.e., 
that i?o[C/]^C^^^^^^^' in the Thomas-Fermi Umit of large 
U. Also, Eq. (O shows that the excitation energy with a 
breathing-mode oscillation of amplitude p scales as with 
the oscillation amplitude p; we will use this for our energy 
analysis below. 

Radius dynamics interpretation — We now proceed to explain 
the Q{t) behaviors presented above, in terms of radius dy- 
namics. Figure |2] (top row) shows the radius evolving as a 
function of time for various ramp shapes, for reasonably large 
T. In the center and bottom rows, we show the deviation of 
a{t) from the equilibrium radius corresponding to the instan- 
taneous value of the interaction, Ro{t) = Ro[U{t)]. For a 
truly adiabatic ramp, a{t) would follow Ro{t) exactly; there- 
fore the deviation f{t) = a{t) — Ro{t) is at the heart of non- 
adiabaticity and the behavior of this quantity determines the 
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FIG. 2: Top row: condensate radius dynamics a{t) for various 
ramp shapes, r — 15. Center and lower rows: deviation f{t) from 
tlie 'instantaneous' ground-state radius Ro(t). Dimensions D — 2, 
{U„Uf) = (20,21). 



amount of heating. After the ramp, the f{t) function is purely 
oscillatory; the heating scales as the square of the oscillation 
magnitude. Figure |2]presents radius dynamics for D = 2; the 
ID and 3D cases are very similar 

For the [A] ramps with derivative discontinuities (middle 
row), the magnitude of the final oscillations of f{t) is deter- 
mined at the ramp endpoint. For the (t/r)" ramp, the os- 
cillation magnitude is 0{t~") during the ramp, and turns 
into 0{t~^) after the endpoint kink. For a > 1, the final 
0{t^^) oscillation is parametrically larger than the during- 
ramp 0{t^°') oscillation. 

We will first explain the ~ r^^ scaling of oscillations ini- 
tiated at the kink. If we neglect the smaller oscillations at 
t < T, the radius a{t)^RQ{t) at the kink t = t has "correct" 
value for U = Uf, i.e. / is negligible. However the deriva- 
tive is nonzero, a' {t)f^R'Q{t)\t=T, which scales as ~ r^^. 
Thus we have the following "initial" conditions at t = t+ 
for subsequent evolution: /(r) = 0, f'{T^) = cqt^^. Us- 
ing f{t > r) w psm{iU + S), these initial values imply 
p ^ . This explains the 0(r~^) oscillation magnitude 
and hence (!)(t^^) residual energy for ramps having a deriva- 
tive jump at the endpoint. 

The oscillations of Q{t) (Figure [TJl) can be explained by re- 
laxing the approximation f{t < t) w made above. The 
small oscillations of f{t < r) guarantee that a'{t = t) os- 
cillates around i?o(< = r) as a function of r. This results 
in the final breathing mode amplitude p to oscillate around 
its 0{t~^) value as a function of r. Since the f{t < t) 
breathing-mode strength is smaller for larger a, the oscilla- 
tions of the heating with r are weaker for larger a, as seen in 



Figure[TJl. 

The lowest-row panels of Figure |2] focus on the smooth [B] 
ramps. In these cases, the breathing-mode strength (~ t~") 
initiated at the beginning of the ramp remains unchanged; 
there is no kink to abruptly create larger oscillations. We 
therefore need only to explain the strength of oscillations at 
the beginning of the ramp, where r{t/T) = 1 — e^''/"^)" w 
(t/r)". 

We first rewrite Eq. (|2]i as an equation for f{t). To simplify 
notation, we will write this out explicitly only in the Thomas- 
Fermi limit, Uij ^ 1, and small oscillations, /(t)<Ci?o(0- 
(The arguments can of course be modified to go beyond the 
Thomas-Fermi restriction. Small f{t) is guaranteed for large 
r.) We obtain 



{D + l)u'{t)^ 

2D + 3 

(L» + 2)2mT^ 



= 0, (4) 



with u = U/{2tt)-^/^. The first two terms give pure oscilla- 
tory behavior (breathing mode at fixed u); the last two terms 
are corrections due to time-varying interaction. 
We first treat ramps with zero initial slope, i.e., a > 1. The 
initial conditions at t = 0+ are then /(O) = /'(O) = 0. With 
u ~ Ui + {5u){t/T)", the u" correction is dominant compared 
to the u'^ correction att t. The dominant correction terms 
take the form ci/r^ for a = 2, and cit/r^ for a — 3. The 
solutions of the resulting differential equation are linear com- 
binations of oscillatory trigonometric functions and algebraic 
functions. It is straighforward to verify that the boundary con- 
ditions /(O) = /'(O) = force the oscillatory part to have 
coefficients scaling as ~ t"^ for a = 2 and ~ t^^ for a = 3. 
This explains the Q ^ r^^" behavior for integer a > 1. 
The a ~ 1 case is slightly different. The initial condition 
still involves ct'(O) = 0, but since Ro{t) = [u(t)]i/(^+i) has 
finite slope at t = 0+, this now corresponds to /'(0+) = 
— i?o(0^) = —cj/t. This initial condition leads to a 
purely oscillatory f{t) with amplitude ~ t^^, which explains 
Q{t) - fora = 1. 

Beyond pure radius dynamics — We justified our single- 
parameter analysis by noting that Q{t) /Q{0) is reproduced 
splendidly by such a description (Figure [T];). The GP dy- 
namics is of course richer than this minimal description, one 
indication of which is that the un-normalized heat function 
Q{t) obtained from Eqs. (O deviates from full-GP re- 
sults. Although a complete study of all aspects of shape dy- 
namics induced in a ramp is beyond the scope of the present 
investigation, in Figure [3] we show some basic additional ef- 
fects. After the size, the next obvious shape characteric is the 
kurtosis k, related to the fourth moment of a distribution, such 
that a gaussian has k = and a more "rectangular" (sharper- 
peaked) distribution has k < (k > 0). The left panels 
of Figure [3] show that in a slow ramp, the kurtosis oscillates 
around the instantaneous equilibrium value just as the radius 
does. In the right panels, the large-amplitude radius oscilla- 
tions after a sudden (r = 0) quench is seen to vary from the 
single-parameter description; there are dramatic jumps of K{t) 
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FIG. 3: Physics beyond pure radius description. D = I, (Ui,Uf) 
= (0,10). Left panels: r = 15 ramp. Right panels: instantaneous 
quench. On top panels, radius from full GP calculations (full lines) 
compared to single-parameter variational calculations (dashed lines). 
Lower panels show kurtosis, not accessible in the radius-only de- 
scription. 



at the times when the (j{t) deviation is prominent. The en- 
ergy clearly leaks from breathing-mode oscillations into other 
channels, even though the total energy is conserved. 
Relation to other results; Open questions — We have ad- 
dressed the fundamental question of adiabaticity in the con- 
text of a paradigm system of cold-atom physics, namely, a 
Bose condensate in a trap. We find that the heat function 
Q(t), which characterizes non-adiabaticity, to have an over- 
all power-law decay. The general intuition is that the decay 
of Q{t) should be exponential if there are no gapless points 
in the path of the ramp; thus finite systems are generically 
expected to have exponential decay, e.g., Ref. |0]. Remark- 
ably, the GP description successfully mimics a gapless ther- 
modynamic limit by providing power-law decay of Q{t), al- 
though concepts like gap or density of states are not meaning- 
ful within the GP description. 



Our result Q{t) 



for smooth ramps is consistent with 



adiabatic perturbation theory: if our ramps are put into the 



form r{t)^vt°' / a\, we get Q{t)^t °'^v , which is the 
generic perturbative expectation f^. In the formulation of 
Ref. 101 (Sec. 3), the F{x) function for our ramp can be shown 
to have asymptotic form which translates to an "extrin- 

sic" contribution Q{t) ^ r^^". The GP description thus re- 
tains nontrivial dynamical information pertaining to the full 
quantum description, despite being "merely" a nonlinear dif- 
ferential equation. 

Another feature we have explored is the sensitivity to a final 
kink in the ramp shape. A recently discovered effect of such 
kinks is logarithmic contributions to Q{t) The 
effect we have found (kink induces larger oscillations over- 



whelming initial excitation) is quite different. It is an open 
question whether or not this is unique to the present system. 
Oscillations of Q{t) are relatively poorly understood, and 
may well be generic in many-body ramps. In our case, it ap- 
pears explicity due to ramp shape kinks. Like other Q{t) fea- 
tures, we have provided a very physical interpretation in terms 
of radius oscillations. In other known examples of Q{t) oscil- 
lations 0-|6l], the physical explanation of oscillations, where 
known, are all different. 

The present work opens up several new research avenues. 
Ramps in the trapping frequency should also induce radius 
oscillations, but details may well be different from interac- 
tion ramps. Physical insights developed in our study of Q{t) 
can perhaps be applied to better understand "optimal ramp" 
studies seeking to find ramp paths producing minimal heat- 
ing 111 111 . Another natural extension of our work is a treat- 
ment of condensates beyond the GP description, e.g., trun- 
cated Wigner schemes or numerical full quantum treatments 
of few-boson systems. 
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